EDGE STRESSES IN A COMPOSITE PLATE WITH A HOLE 
- A HYBRID FEM APPROACH 


A Thesis Submitted 

in Partial Fulfilment of the Requirements 
for the Degree of 

MASTER OF TECHNOLOGY 


by 

MAKARAND JOSHI 


CO the 

DEPARTMENT OF MECHANICAL ENGINEERING 

INDIAN INSTITUTE OF TECHNOLOGY KANPUR 


JULY. 1*91 



f 2 NOV '«>91 
(661 ' ■' f 

central library 



4oe.Jto. 



CERTIFICATE 


O; 

i O 

ft. 1,1. r K;.ti^ 
\ Siibmi(T.C'i "(1.1^1- 

V. 




This is to certify that the work entitled "Edge Stresses in < 
Composite Plate With a Hole - A Singular Hybrid FEM Approach" by Hakaranc 
Ganesh Joshi has been carried out under my supervision and has not beer 
submitted elsewhere for a degree. 


i'^«r 


N. N. KISHORE 




Assistant Professor 

Department of Mechanical Engineering 

Indian Institute of Technology 


KanpLff"; 208016 



1 


ACKNOWLDGEMENTS 


I am extremely grateful to Dr. N.N. Kishore for his able guidance, 
critical appraisals and thought provoking suggestions, from the inception 
to the completion of this work. 

I feel indebted to Vijay Pankhawala, Vinu and Hadhav who helped me 
draw the figures and graphs. I would like to thank all in the ESA lab. who 
helped me take the rough and the final prints of the thesis. I would also 
like to thank GX. Singh and Subrato Saha for the enlightening discussions 
on this work. A million thanks to the ghat 'Aish'ociation, GK, Shankh, Rajiv, 
Nachiketa, Kak and all others who made my stay at UT-K a pleasant and 
memorable one. 


Makarand Joshi 



C04TENTS 


PAGE 

CERTIFICATE .. i 

ACIO>«WLEDGEMEhfTS .. n 

CONTENTS iii 

LIST OF SYhffiOLS .. v 

LIST OF RGURES .. vii 

SYNOPSIS .. X 

CHAPTER I INTRODUCTION 

1.1 Introduction .. l 

1.2 Literature Survey 2 

1.3 Scope of present work .. 6 

CHAPTER n PROBLEM DEFINITION 

2.1 Introduction .. 8 

22 Problem definition .. 8 

CHAPTER m TEE SINGULAR HYBRID FEM APPROACH APPLIED TO 
LAMINATED PLATES 

3.1 Introduction .. 12 

3.2 Hybrid stress elements .. 13 

3.3 Layer stress interpolation for non-singular 

hybrid and singular hybrid elements .. 22 

CHAPTER IV RESULTS AND DISCUSSIONS 

4.1 Interlaminar stresses near the circular hole 

boundary in symmetric composite laminate .. 39 

4.2 Interlaminar stresses near the square hole 

boundary in symmetric composite laminate .. 42 



CHAPTER V CONCLUSIONS AND SCOPE FOR FUTURE WORK 


i'w 



5.i 

Conclusions 

.. 67 


5.2 

Scope for future work 

67 

REFERENCES 



.. 69 

APPENDIX A 



.. 73 

APPENDDC B 



.. 76 



LIST OF SYMBOLS 


X^, X2, X3 
X, y. z 

e 


Material co-ordinate system 

Global co-ordinate system 

Angle between the material and global 


t 


Vni 



I 


^11' ^22 


^ 12 ' ^23 

Nj(f . n) 


co-ordinate systems 
'v^ector of stresses in layer i 
Vector of strains in layer i 
Displacement vector for layer i 
Compliance matrix 
Prescribed tractions 
Volume of layer i 

Portion of element boundary on which tractions 
are prescribed 

Total number of layers in the laminated 
composite 

Young's moduli in the fibre and transverse 
directions 

Major and minor Poisson's ratio respectively 
Matrix of polynomials of stress assumption 
Constants used in stress assumptions 
Biquadratic Serendipity shape functions 
The normalised co-ordinate system 
Nodal displacements of layer i 
Shape functions in normalised co-ordinates 
Strain displacement differential operator matrix 



I J I 

ti 

K 


Q 

Cs^ 



X 


r, e, z 


s,^ 

•"eg 


f'en 


OTp 


o-z 

*^0 


'^rz 

H 


Jacobian for co-ordinate transformation 
Half the thickness of layer i 
Qemental stiffness matrix 
Qemental load vector 

Matrix relating layer stress parameters to 
laminate stress parameters jS 
Matrix relating nodal displacements to 
laminate 

nodal displacements q 

Power of singularity 

cylindrical co-ordinates 

polar co-ordinate radius and angle 

Radius of centroid of the sector 

Radius of the mid-node on the boundary of the 

sector 

Stress in the radial direction 

Stress in the axial direction 

Stress in the circumferential direction (Hoop 

stress) 

Shear stress 

Strain in the radial direction 
Strain in the axial direction 
Circ^umferential strain 
Shear strain 



Vll 


LIST OF FIGURES 

RGURE NO. TITUE PAGE 

2.1 Geometry of laminate with a circular hole .. 10 

2.2 Geometry of laminate with a square hole .. 11 

3.1 Geometry and layer numbering of multilayer 

plate element .. 33 


32 Lamina with its principal material axes 

oriented at angle O with reference to the 


global axes .. 34 

3.3 Displacement degrees of freedom .. 35 

3.4 Singular elements at hole boundary .. 36 

3.5 Mapping from natural coordinates S/ ^ to 

physical coordinates s and 0 ..37 


3.6 Position of any point on sector w.r.t. axisymmetric axes. 38 

4.1 FEM Discretization of laminate plate with circular hole 


<85 eleiments) .. 45 

42 FEM Discretization of larriinate plate with circular hole 

(17D elements) .. 46 

4.3 Angular variation of power of singularity .. 47 

4.4 cTz radial distribution for 0/80 laminate with circular 

hole (85 elements) e = 0° ..48 

4.5 tj2. radial distribution for 0/90 laminate with circular 

hole (85 elements) 0 = 45° .. 49 

4.6 cTz radial distribution for 0/90 laminate with circular 

bole (85 elements) 6 = 90° .. 50 

4.7 radial distribution for 0/90 laminate with circular 



Vlli 


hole (170 elements) 6 = 0° ..51 

4.8 ( 7 z radial distribution for 0/90 laminate with circular 

hole <170 elements) 6 = 45° ..52 

4.9 Cz. radial distribution for 0/90 laminate with circular 

hole <170 elements) 6 = 90° ..53 

4.10 Tpz radial distribution for 0/90 laminate with circular 

hole (170 elements) .. 54 

4.11 00 through thickness distribution for 0/90 laminate 

with circular hole <170 elements) .. 55 

4.12 Oz thirough thickness distribution for 0/90 laminate 

with circular hole (170 elements) .. 56 

4.13 Oz radial distribution for 90/0 laminate with circular 

hole <85 elements) .. 57 

4.14 Tpz radial distribution for 90/0 laminate with circular 

hole <85 elements) ..58 

4.15 (70 through thickness distribution for 90/0 laminate 

with circular iTole (85 elements) .. 59 

4.16 (7z through thickness distribution for 90/0 laminate 

with circular hole (85 elements) .. 60 

4.17 FEM Discretization of laminate plate with square hole 

<85 elements) .. 61 

4.18 FEM Discretization of laminate plate with square hole 

<170 elements) .. 62 

4.19 (7z radial distribution for 0/90 laminate with square 

hole <85 elements) .. 63 

4.20 Tpz radial distribution for 0/90 laminate with square 


hole <170 elements) 


64 



iX 


421 


4.22 


through thickness distribution for 0/90 laminate 
with circular hole (85 elements) 

CTz through thickness distribution for 0/90 laminate 
with circular hole (85 elements) 


.. 65 


.. 66 



SYNOP^ 


Composite materials have found wide application in the aerospace 
and aircraft industry, due to their decisive advantages in terms of 
structural properties, over the other conventional materials. This has 
boosted the research on such materials in the past few years, mainly in 
the areas related to the manufacture, and fracture and fatigue. 

Due to the design requirements, the composite laminates may have 
holes, cut-outs etc. These provide free edges where, a three dimensional 
state of stress develops even under uniaxial loading. This happens due to 
the sudden change in the material property at the interface. Delamination 
is a major cause of failures in composite laminates. The present work, 
endevours to analyse the edge stresses in symmetric crossply laminates 
having central circular/square hole, at the hole boundary, which is a 
potential region for the delamination ttfloe of failures. 

Many attempts have been made in the past to solve this problem by 
the application of semi-analytical and finite element methods. The 
analytical methods are complex in their formulation and solutions, while the 
displacBiTvent based FEM solutions require a very fine mesh and large 
computing times. Hybrid FEM solutions, which have been developed in the 
recent past, provide expeditous computation but do not guarantee 
convergence at the edges of interfaces. The use of singular elements can 
overcome this drawback. The appropriate singular elements, when used in 
the regions close to the hole model the singularity of the stresses and 
improve convergence. 

In the present case, a l%4>rid stress 8-node isoparametric multi- 



layered plate element which includes an appropriate singular term has been 
used. The hybrid stress model is based on the modified complimentary 
energy principle and takes into account the transverse deformation 
effects. All the displacements are assumed to vary linearly through the 
thickness. 

The singular element^ which is used in the boudary layer region uses 
the hybrid stress interpolation in conjunction with a singular term. The 
power of singularity has been found by solving a separate eigen value 
problem. 

A computer code in FORTRAN language has been developed to solve 
the above problem. The results show significant improvement in the 
accuracy and efficiency of the present method compared with the earlier 
methods. Also, they indicate a mixed mode failure near the hole boundary. 



CHAPTER 1 


imRfjDucrnaN 


1.1 Cofnposite Materials : 

Composite materials are made of two or more constituents; which on 
a macroscale have distinct interface separating them. These materials 
achieve certain physical properties not realizable in their constituents; 
such as high specific strength and stiffness; and controlled anisotropy. As 
they are usually endowed with good fatigue and corrosion resistance; low 
heat conductivity; good electrical insulation properties, favourable cost 
effectiveness etc.; they find wide applications in aerospace and transport 
industry. In recent past extensive research work has been carried out on 
the manufacture; characterization; fatigue and fracture of these materials; 
to promote further use of these materials. 

1.1.1 Interlaminar Str^ses/Edge Effects : 

A prominent disadvantage of laminated composites is that under 
simple uniaxial or biaxial loading, they develop a triaxial state of stress 
near the free edges. These out-of -plane stresses are called the '■'inter- 
laminar'^ stresses and this phenomena is called the edge effect, as it is 
evident only in the close vicinity of the free edge, in a region 
approximately equal to the thickness of the laminate. These interlaminar 
stresses are of critical concern to designers as they may lead to 
delamination type of failures at loads well below those corresponding to 
inplane failure. Interlaminar stresses may be present as a result of 
mechanical; thermal or moisture loading. 
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These stresses between layers are caused due to the mismatch of 
properties of bonded adjacent layers and the presence of gradients in the 
inplane stresses. Stress free boundaries present the most severe stress 
gradients and hence the most severe interlaminar stresses. 

In contrast to the isotropic materials, the anisotropic materials 
present a very broad range of variation in properties across the adjacent 
layers. The material properties which result in a mismatch of properties 
are primarily the Poisson's ratio, the coefficients of mutual influence, and 
the thermal and hygroscopic expansion coefficients .The transverse shear 
stresses result in the sliding of adjacent layers over one another. Though 
this phenomena is localised near the free edges and insignificant in the 
far field regions, it may prove fatal under fatigue loading as it may lead 
to matrix cracks. Also, the transverse normal stress may lead to 
delamination type of failures. 

E>ue to functional or design requirerrents, the composite may contain 
htoles, cut-outs etc. These, apart from the free edges show stress 
singularity. Hence the analysis of laminated composites to estimate the 
interlaminar stresses proves important from the design point of view. 

1.2 Literatix^ Survey : 

The existence of interlaminar stresses along the free edges, in the 
boLTKiary layer region of a laminated composites has been known for some 
time. The first pL±)lication dealing with the singular behaviour of stresses 
at the intersection of a boundary and bonded dissimilar materials is that 
of Bogy (1968). The earliest work dealing with the anisotropic materials 


appears to be the study by Puppo and Evenson (1970), in which, the laminate 
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was modeled as a set of anisotropic layers separated by an isotropic layer 
of adhesive material having negligible stiffness. 

The first complete three-dimensional analysis of interlaminar 
stresses by applying finite difference technique was done by Pipes and 
Pagano <1970). They showed that interlaminar stresses existed along the 
free edge in the boundary layer region and that its width was 
approximately equal to the laminate thickness. The singular behaviour of 
the transverse normal stress Oz and the transverse shear stress Tpz 
varies to a large extent with laminate configuration (material, fiber 
orientation, laye^r thickness and stacking sequence). These results were 
later experimentally verified by Pipes and Daniels (1971). 

The relationship between interlaminar stresses and delamination was 
analytically disscussed in many subsequent papers. Pagano (1974) presented 
an approach which accounted for the influence of pertinent material and 
geometric parameters on the distribution of normal stress cr^ along the 
central plane of a symmetric finite width composite laminate. Pipes and 
Pagano (1974) developed an approximate elasticity solution for the 
response of finite width, symmetric angle-ply laminate mder uniform axial 
strain. 

Tang (1976) derived the governing differential equations for 
uniformly loaded rectangular composite plate using the classical laminate 
theory and approximate boundary conditions. Tang (1977) later used a 
boundary layer theory to obtain stress distribution around a circular hole 
in C0/901s and [±45]^ laminates, using very large values for the ratio of 
hole radius to ply thickness. 


Pagano (1978) reduced the free edge class of boundary value problem 
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to a one-dimensional oroblem to predict the stress fields. The solution 
involved consideration of 13N algebraic and ordinary differential 
equations, where N stands for the number of layers. 

Kassapoglous and Lagace (1986) presented a method which predicted 
interlaminar stress state. The expressions for stresses were in terms of 
exponentials based on the shape that the interlaminar stress must take in 
order to ensure overall force and moment equilibrium. Boundary conditions 
and traction continuity between plies were satisfied exactly. 

Chaudhari and Siede (1987) presented an approximate semianalytical 
method to determine transverse shear stress variation across the 
thickness. The method utilised a quadratic displacement potential energy 
based FEM. 

A three-dimensional finite differnce model by Altus et^I. (1980) 
followed a finite element formulation by Isakson and Levy (i^i) to 
determine the interlaminar stresses for generalised plane strain 
conditions based on the work of Pipes and Pagano (1970). Subsequently many 
results were published using finite element techniques. 

Rybicki (1971) presented a three-dimensional FE formulation using 
complementary energy approach. Wang and Crossman (1977) included 
mechanical, thermal and hygrothermal loads, in their FE formulation for a 
symmetric finite width laminate under tniaxial tension. This formulation 
was based on the generalised plane strain approach. Wang and Dickinson 
(1976) applied the extended Galerkin method with displacements and 
interlaminar stresses satisfying the geometric boundary conditions 
represented by a complete set of Legendre polynomials. The final solutions 
were obtained by requiring the satisfaction of continuity conditions at 
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each interface and the stress boundary conditions at the exterior planes. 

Lee (1980) adopted three de^ees-of-freedom per node for three-dimensional 

FE analysis of symmetric composite laminate with a center hole and 

subjected to biaxial inplane loads. Raju and Crews (1982) investigated edge 

stresses near the boundary of a circular hole in a symmetric cross-ply 

laminate using a quasi-three-dimensional FE analysis with 19000 degrees- 

of -freedom, i.e., a very fine mesh was used in the boundary layer region 

which provided more details in the neighbourhood of the singularity. They 

confirmed the presence of stress singularities and determined the power 

* 

of singularity and the strength of singularity by assuming 

C7=:Ar‘^ 

where A stands for the strength of singularity and \ the power of 
singularity. They found these to be same for a number of ply combinations. 

Pian (1973) developed the hybrid element and showed that because of 
the high accuracy of the h^orid elements, they are more efficient than the 
assumed displacement elements. 

In the recent past, Spilketr et.al. have developed stress formulation 
to analyse various types of plate problems based on different assumptions. 
Spilker and Chow (1980) analysed symmetric cross-ply laminates under 
uniform in-plane strain using the hybrid approach, in which, equilibrium in 
each layer, continuity of appropriate stresses and displacements at 
interface between layers and the traction free boundary conditions were 
satisfied exactly. They disscussed the importance of the satisfaction of 
the traction-free edge conditions. Non-satisfaction of these conditions 
could lead to a spurious or oscillating result in the boundary layer 


region. Spilker and Munir (1980) developed an assumed stress FE model for 
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the analysis of thin plates. They carried out a study for determining the 
criteria for the selection of stress assttfnptions to avoid element locking 
in the thin plate limit. They later (Spilker and Munir 1980) developed a 
hybrid stress element for the analysis of thin and moderately thick plates. 
They showed that the hybrid stress element yields better accuracy as 
compared to the assumed displacement based Mindlin plate element with 
reduced integration. Spilker (i960) developed a hybrid element for thick 
laminates wherein attention is restricted to the cylindrical bending . 
Spilker (1982) also developed a hybrid stress 8-node isoparametric 
multilayer plate element for the analysis of thin to thick fiber reinforced 
composite laminates. 

Analytical solutions obtained by the use of complex potentials and 
eigen function series by Wang and Choi (1982) were used by Wang and Yuan 
(1983) to develop a singular hybrid element . This element, in conjunction 
with the conventional isoparametric quasi-three-dimensional elements has 
shown to be suitable for examining the detailed nature of boundary layer 
stresses in composite laminates. 

1.3 Scope of Present Mark : 

The purpose of the present work is to analyse the edge effects in a 
symmetric , finite, multilayered rectangular composite laminate having a 
central circular/square hole. The present work uses the hybrid element 
developed by Spilker (1982) in conjunction with the singular el&nent 
developed by Suresh (1987). The singular element is used in the boundary 
layer region while the hybrid element is used in the far field region. The 


exact order of singularity has been found by solving an eigen value 
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problem as done bg Yamada (i9B3) and Suresh (1987). This has been used to 
construct singular hgbrid elements while simpler hybrid stress variation is 
used for the far field elements. 

The aim of the work is to ensure a computationally efficient method , 


by using a coarse mesh. 



CHAPTER 2 


PROBLEM DEFMTION 


2.1 Introduction 

In the present work, a composite plate with a hole has been analysed. 
The problem deals with the edge stresses in [0/903s and ESO/Ols laminate 
having a central circular /square hole. Due the complex nature of stresses 
a hybrid FEM approach has been employed. To simulate the edge 
singularities, a -singular element approach has been used in the region 
close to the boundary. 

22 . Problem Definition 

Figure2.i shows a four ply laminate made of graphite/epoxy with 
each ply having thickness h. It has a central circular hole of radius Pq. 
Figure2.2 shows a similar laminate with a central square hole. These 
laminates are subjected to an uniform far-field displacement in the y- 
direction as shown. Each ply is idealised as homogenous, elastic, 
orthotropic material having following properties -■ 

= 20 X iO^ Psi (138 GPa) 

E 22 = 2.1 X 10 ^ Psi (14.5 GPa) 

G^ = 0.85 X 10^ Psi (5.86 GPa) 

G 23 = 0.85 X 10® Psi (5.86 GPa) 

“ 0.2i 


^23 ~ 0 . 2 i 
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Here, the subscripts 1,2 and 3 refer to the material axes directions. 
Axis 1 is in the direction of the fibers and axis 2 is in the plane of the 
ply perpendicular to axis 1. Axis 3 is in the direction normal to the plane 
of the laminate. 

The analysis has been carried out for two layup configurations viz. 
C0/90 ]e and [90/015. Because of the symmetry in the problem, only one-eight 
of the laminate was analysed with the following boundary conditions : 

.i ) Displacement boundary condition 

a) u = 0onx = 0 

V = 0 on y = 0 

w = 0 on z = 0 

b ) Dn the far end y = 1 plane 
u = uo ( constant ) 
ii ) Traction boundary conditions : 


a) Edges x = b, the top surface 2 = 2h , and 
the hole boundary r = rO 


Traction Free 


The analysis is aimed at obtaining the stress distribution in the 
idealised laminate caused by the above mentioned loading with emphasis on 
the edge stress near the boundary. In the present work singular elements 
in conjunction with the hybrid finite element method has been used to 


obtain efficient solutions. 





CHAPTER 3 


THE SB^WULAR HYBRID FEM APPROACH APPLIED TO LAMINATED PLATES 

3.1 Introdijction : 

In problems of linear elasticity involving an interface joining two 
dissimilar materials, often complex states of stress are encountered. This 
is due to the existence of stress singularities in a localised region close 
to the free surface and at the interface. To obtain satisfactory results of 
such problems, as in the present case, a hybrid FEM approach in 
conjunction with the application of singular element in the boundary layer 
is the appropriate method. 

In linear solid mechanics, the principle of minimum potential energy 
involves a compatible displacement field as its variable with the stress 
equilibrium equations and the traction boundary conditions as its Euler - 
Lagrange equations. The complimentary energy principle involves an 
equilibrated stress field satisfying the traction boundary conditions, as a 
variable with equations of kinematic ccmpatability and displacement 
boundary conditions as its Euler - Lagrange equations. 

For the solution of a linear boundary value problem by the finite 
element method by the application of the minimum potential energy 
principle, the displacement field should not only be continuous within the 
element but also be compatible at the interel^ment boundaries. The 
variational principle then gives (a) equilibrium within the element and (b) 
the traction reciprocity condition as the Euler-Lagrange equations. 
Similarly, for the application of complimentary energy principle, the stress 


field should satisfy the equilibrium equations and the interelement 
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traction reciprocity conditions. The Euler-Lagrange eouations in this case 
are (a) the kinematic compatability within the element and (b) interelement 
displacement compatability. 

For problems related to plate bending, shells, multilayered 
composites, problems with singularities, etc the application of the primal 
principles to FEM restrict the flexibility and do not give satisfactory 
solutions. To obtain this flexibility, one may relax the conditions of 
displacement continuity and/or traction reciprocity at the interelement 
boundaries on the admissible displacement/stress fields and introduce 
these as posteriori constraints by use of the Lagrange multiplier 
technique. This is the underlying concept of the Hybrid FEM methods. 

If conventional displacement methods are used for the solution of 
problems involving stress singularities, a large number of elements have to 
be used to obtain credible results. On the contrary, the application of 
hybrid element formulations can ensure better results. However, even these 
methods do not guarantee convergence. This can be overcome by the use of 
appropriate singular element in conjunction with the hybrid FEM 
formulation, which will ensure efficiency in computation, and convergency 
of results. 

32 Hybrid Stress Element 

These elements are refered to as the non-singular elements in the 
ensuing discussion. The elements considered here are assumed to lie in the 
x-y plane with the z-axis corresponding to the transverse direction. The 
element is composed of I perfectly bonded surfaces Qnterlayer surfaces 
and lower/upper surfaces) which are numbered from the bottom to top of 
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the laminate. The origin of the reference surface, z=0, is at ttre bottom and 
the location of each surface is given by 2 =hj^, (h^*0), h 2 > ■■■■i hj^j^ <see fig 
3.1.). Each layer is of constant thickness and constant material properties 
(i.e. material constants in the material axis system, and ply/fiber 
orientation angle, 6), but thickness and properties may vary from layer to 
layer. 

The elements used in the present problem are based on the hybrid 
stress model. The hybrid stress model uses a modified complimentary 
energy principle in which the interelement traction continuity and the 
mechanical boundary conditions have been relaxed via the Lagrange 
multilplier technique. The hybrid stress functional for multilayer element 
can be expressed as (Spliker 1982) 

= SZ { 4 / 2 ^^ i" dV . / 2 ^^ dV + / f dS } (3.1) 

^ ^ ^ 

where 2 ^ is the vector of stresses for layer i, is the vector of strains 
in the layer i computed (via strain-displacement relations) from 
displacements, u^, is the Compliance matrix of the material ( 2^ ^ 

where are actual strains, to be distinguished from ), are 
prescribed tractions, is the volume of the layer i for element n, and 
ScTn is the portion of the element boundary on which tractions are 
prescribed. The inner sum is over the layer I, while the outer sloti is over 
the number of elements. 

Each layer of the multilayer element is an tanidirectional lamina and 
for the present purpose, is assumed to be homogenous and orthotropic 
Material properties are specified with respect to the material axes (i.e. 
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co-ordinates X2< Xg where Xj^ corresponds to the fiber direction, X2 is 
the in-plane co-ordinate perpendicular to Xj^, and Xg is the transverse co- 
ordinate, refer fig. 3 . 23 . The required compliance matrix relates 

stresses and strains in the global co-ordinate system (x,y,z). The material 

and global co-ordinates are related by the angle 6, where 6 is the angle 

measured from the global axis x to the fiber direction X^^. Applying 

transformation laws and assuming stresses in the i^^ layer as 
:T 


S^= 


X 

b 

(T^y 0- 

_i 

z yz 

<^^xz 

(^^xy 

1 

Sll 

S12 

®i 3 

0 

0 

Si6 
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®22 

®23 

0 

0 

®26 

®31 

®32 

®33 
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0 

®36 
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0 

0 

S44 

S45 

0 

0 

0 

0 

S54 

S55 

0 


^62 

®63 

S 64 

^65 
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y^ere, with m=cos6 and n=sin6 


c — 

Sil - r 


ii 


E. 
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®13 = 
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E 


11 
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22 
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1 
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^11 
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In equation (3.2), ^11 and E 22 Yolb^s mcduli in the fiber and 
transverse directions respectively, and 1/23 major and minor 

Poisson's ratios. 

Returning now to the element formulation based on equation (3.1), 
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stresses in the layer are assumed in terms of a finite set of stress 
parameters, for layer i, in the form 

(3.3) 

where - consists of polynomial terms and is defined so that the 
homogeneous equilibrium equations are exactly satisfied. The typical set of 
polynomials used for the present analysis are given in Appendix-A. 

The isoparametric transformation used in the present 8-noded 
element is 

8 

X = N, (S,r?) X,- (3.4) 

3=1 ^ ^ 

8 

y = S Nj (C.J?) yj (3.5) 

where Nj (J,??) corresponds to the biquadratic serendipity shape 
functions (given in the Appendix-B) and (Xj,yj) are the global in-plane co- 
ordinates of the node j. In addition, a normalised transverse co-ordinate, t, 
is introduced for each layer which is related to z by 

z = I [(h- + h^^^) + T (h^.^.i + h^)3 (3.6) 

so that T takes the value -1 at the bottom of the layer i and +1 at the top 
of the layer i. Using Eqs. (3.4 - 3.6), in the x,y,z co-ordinates can be 

expressed in terms of normalised co-ordinates (J,)?,t). 

The displacements for the layer i are interpolated in terms of 
nodal displacements, q^, for layer i in the form of 

y^ = q^ (3.7) 

where are the shape functions in the normalised co-ordinates (5,7?,t>. 
Here, y^ is to be ordered as y^ = t u v w 3^ (displacements in the x,y and z 
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directions respectively). By applying the linear strain-displacement 
relations to Eq. (3.7), the strains can be related to as 

q^ = q^ « I ^ I B*^ q^ (3.8) 

where is the strain-displacement differential operator matrix, and | J 1 
is the Jacobian of the transformation given by the Eq (3.4 - 3.5). If the 
components of are ordered as in is given by 



D 


0 


( ^ J. 
I ^ 3t? 


3x a 

§17 §S 


0 



0 


0 

1 J I a. 
IT 9? 


I J I a. 

f: SS 


J I a 

i 

0 


q 3? 


f is J. 3x 1. 1 
1 a$ ar? ar? a? 7 


( 8y a 
\ §r? a? 


ay J- ] 
a? ar? ' 


€3 


{ ix ± ix d ) ( ^ a 9yj.i 

v as ar? ar? §s M ar? §5 ■ as ar? ; 


where the transformation of Eq. (3.4 - 3.5) is used in Eq. (3.9) and tj^ is half 
the thickness of layer i, i.e. 



Substituting eqs. (3.3), (3.7), and (3.8) into Eq (3.1) and using the following 

layer matrices •• 

_ 1 1 1 .T . . 

= t^ J J J 1 J I dS dT? dT 

-1 -1 -1 


(3.11) 
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df dr? dr 

-1 -1 -1 


(3.12) 




ds 


(3.13) 


TTmc becomes 

•^mc “ S S { I - )3^ ) (3.14) 

““ - - 7 

At this stage, stresses and displacements are independant from 
layer to layer. However, enforcement of traction continuity aJLong inter- 
layer surfaces will require that certain of the be related to certain 

of the 0^ (i.c, at the interface between layers i and i+i). Similarly certain 
of the nodal displacements, will be related to certain in order to 

satisfy displacement continuity along the interlayer surface between 
layers i and i+1. 


Corresponding to these interlayer surface continiuity conditions, 
layer stress parameters, 0^, and nodal displacements, q^, can be related to 
the overall laminate stress parameters, 0, and nodal displacements q, 
respectively, in the form : 

^ = Cs^ ^ (3.15) 

q (3.16) 

By substituting above equations into Eq.(3.14), 'jcnic 

simplified in terms of laminate/element stress and displacement parameters 
and from Eq.O.l), is given by 

’^mc = E { I H ^ - g'*' G q + q’’’ Q } (3.17) 


where 
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H = 





Cs^ 


(3.18) 


G 



G^cJ- 


G.19) 



(3.20) 


In th. above .Quauon, matrix multiou^ii.,^ are not carried out 
exolicitJu. The operation, in the above epuations are einimilar to element 
aesemblu ooerationsi a set of layer to laminate parameter 'oointer.- and 
nodal disblaoement 'oointery can be used to locate laser matrix 

contributions in the laminate matrices. 


The variation of functional 7rn>c with respect to 5 yeilds the set of 
equations 


H0 - Gq = 0 

Hence, 

0 G q 

and taking variation with respect to q 

g"^ 0 = Q 

Substituting eq.(3.22) into eq.(3.23) 

G^ G q = Q 


(3.21) 


(3.22) 


(3.23) 


(324) 


Rewriting the equation in the form 
k q = Q 

The elemental stiffness matrix can be identified as 
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K = g"’’ G 

*nd Q as the elemental load vector. 


(3.25) 


Once the nodal displacements have been calculated^ stresses •for 
each layer are calculated by obtaining 0 from Eq. (322), j8^ from Eq. (3.15) 
(by extracting the appropriate ^ from §>, and substituting into Eq. (3.3). 
Strains can be obtained from stresses by using S^. 


32.1 Displacement biierpolatian fix' a Layer : 

Moderately thin and thick laminates are characterised by severe 
crosssectional warping; and each layer must be allowed to Lmtfergo non- 

t 

normal crosssection rotations, while preserving interlayer surface 
displacement continuity. To incorporate this warping behaviour, all 
displacement components Cu,v and w in x, y and z directions, respectively) 
are assumed to vary linearly, through the thickness of each layer. The 
degrees of freedom at any node j, which guarantees interlayer surface 
displacement continuity, are therefore u^j, v^^ and where, i signifies the 
points at the bottom, each interlayer and moper surface for this node (i = 

1,2,3 , I+l) (Fig 3.3). The displacement interpolations in normalised 

coordinates for typical la^r i are then given as 

8 

u^ (S,r?,T) = E 2 + 1 

0=1 

8 

(S,T7,t) = 2 Nj (S,T?) c i (i-T) v> + I (1+T) v/'^^ ] (327) 

0=1 

8 

(S,r?,T) = E Nj (S.T}> C I (i-T) W> + I (1+T) ] 

0=1 


( 328 ) 
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where Nj (iJ)) ere the standard S-node biouac^atic Serendipity shape 
functions (shown in Appendix-B). 

The Uj , Sr'S treated as the layer 

degrees-of-freedom, q^ Eqs. (3.26 - 328) then define of Eq. (3.7), from 
which B of Eq. (3.8) can be computed. The ordering of laminate degrees of 
freedom is given by, 


1 .. 2 


^v 2 


1+i. . i. . 2 


q =[Uj^ Uj^ Vj^ 


■••■'^2'*^^ 






the Boolean matrix of Eq. (3.16) can be defined. For 1 layers, the 8-node 
element based on thie 'generar displacement field of Eqs. (3.26 - 328) will 
have 24(1+1) total degrees of freedom per element and 3(1+1) degrees of 
freedom per node. 


3.3 Layer Stress Interpolation for Non-singular Hybrid and Singular Ht^id 
Elements : 

In the present formulation, two types of elements are used. Adjacent 
to the hole boLndary, the stresses are singular and hence singular hybrid 
elements are employed. In the region away from the hole boundary, the non- 
singular hybrid elements have been used. 

3.3.1 Stress Interpolation for Non-singular Elements : 

All components of stresses are included and are interpolated within 
each layer in terms of stress parameters, for that layer, such that 
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homogenous equilibrium equations are exactly satisfied. Introckjcing the 
normalised transverse layer co-ordinatej t, these equations are 


So- 


8x 3y T. 3 t ° 


(3.30) 


8x 8y I. ~Tr = ^ 

i 


(3.31) 


^xz^ acz^ _ - 

9x 8y 8 t 


(3.32) 


The stresses ax^(x.y,T). C7y^(x,y,T), 0 'z^(x,y,T) are asstmed to vary 
linearly in t within the layer. Then from Eqs. (3.30) and (3.31), o-xz^(x,y,T) 
and o-yz^(x,y,T) must be of the order- from Eq. (20.3) (rz^Cx.y.r) must be of 

3 

the order t within the layer. A 67 term (67 6’s), generalised hybrid stress 
assumption used in the present analysis is presented in Appendix-A. It may 
be noted that the stretching portions of o'y^ <^xy^j remain complete 

cubic polynomials in x, y, o-^z^ and o-yz^ include all quadratic x, y terms 
and CTz^ is linear in x and y. 

The layer stress fields satisfy the homogenous equilibrium 
equations within the layer, but stress parameters are independent from 
layer to layer. In the the present element used, the interlayer traction 
continuity are satisfied by establishing a set of element/laminate strsss 
parameters, $, and relating to ^ as in Eq. (3.15). 

The interlayer surface traction continuity requires that O’xzi t^yz 
and Oz t)e continuous at interlayer surfaces. For example for the 
polynomials assumed in the present analysis (Appendix-A), evaluation of 
these stresses at the lower surface (t = -1) of layer i yields : 
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o-^xzIt^-I “ + %^3cy + Pg^y2 ) 

c^^yzlra-i = h ( 07^ + 00^x + 09^y + + ^^^^xy + 0i2^y^ ) (3.32) 

‘^^zl-rs-i = < 0i3^ + 0j,4^x + 0j^5^y ) 

and the upper surface (t = i) of layer i : 

xzIts+i = ^ 053 + %4 ^ '*' %5^y ■*■ ■*■ "*■ ^ 

^ yz^Tss+i ~ *'%9 ■*■ ^60^* ■*■ ■'■ %2^^^‘'‘ %4^y^ ^ (3.33) 

^^zt-j-s+l = < 065^ + ^66^^ ■*" %7^y ^ 


With the simple form of transverse shear and rcrmal stresses at 
interlaye;^ surfaces, the traction continuity condition at the interface 
between layers i-i and i (where i > i) requires that 


a 1 _ ^i-1 a i-i 

- X“ ® j+52 
^i 


j = 1,2, ,12 


(3.34) 


. t2 

a 1 _ i-1 a 

^0 - 71~ ^0+52 
^ i 


3 = 13,14,15 


(3.35) 


thus the set of independent laminate stress parameters, $, for the present 


element is reduced to 
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^^67 ^ 16^ 17 ^^67 *^^6^17 ^67 ^ 

Thus the total number of indspendant stress parameters for I layer 
laminate is (67 - 15) I. 

The Cg of Eq (3,15) can be constructed for each layer using Eqs. 
(3.34 - 3.35) and the laminate 0 of Eq (3.36). In practice, matrix 

multiplications in Eq. (3.17) and (3.18) involving are not performed. 
Instead, row and column and rows in corresponding to sire 

scaled as per Eqs. <3.34 - 3.35), after which a set of layer-to-laminate 
stress parameter ^pointers^ are used to appropriately assemble ttffi scaled 
into H and G into G^. 

3.3,2 Stress biterpolation for the Singular Elenieni : 

In the present problem of analysing crosiRDly comcxDsite laminates 
with a central hole, subjected to far-field loads, the stresses near the 
boundary will be very large and three-dimensional in nature (theoretically, 
some of the components are singular). As the singularity is due to the 
interaction of the free edge (circular hole bomdary) and tt^ material 
discontinuity between layers, it can be observed that the severity of the 
stress singularity changes from point to point along the hole boundary. 
Such a problem is truely three-dimensional in nature and no work has been 
reported so far to solve this problem. However, work was reported 
(Yamada(i983) and SureshdSS?)) wherein, the severity or the power of 
singularuity was determined for an axisymmetric cases. 

To determine the singularity as a first step and its use in the 
hybrid formulation, a simplified procedir'e along with a few assumptions is 


adopted. 
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- At each poimt along the hole boundary; the radial and tangential 
(r;6) material properties are determined and the singularity index at that 
point is found under the assumption that the same conditions prevail all 
along the boundary - thus leading to the axisymmetric nature. As can be 
observed; the following two factors are neglected viz. 

i- 'I'be variation of the material properties along the tangential direction 
near the hole boundary; and 

2. The non-axisymmetric nature of the farfield loads. 

The effect of these two factors in the stress singularity may be assumed 
to be marginal because of the fact that stress singularities arise mainly 
due to the free edges (intersection of the hole boundary and the material 
discontinuity in the thickness direction) and the axisymmetric assumption is 
used in regions close to the hole boundary. Also, the power of singularity 
can be very easily determired for an axisymmetric case. This stress 
approximation is then required to satisfy 


ScTp B7rz , _ 

ar 8z r 

aTrz Sq’t '7'rz 

ar az ^ 


(3.37) 

(3.38) 


As the interlaminar stresses Tpz and the in-plane stress show 
singularities, the following approximations have been used for these 


stresses and the other stress components are then found by enforcing the 


stress equilibrium conditions. 


0-2 = ( ( r - To ^ (3-39) 

Trz = ^3 ( C r - To ^ (3.40) 

(Tq = ( ( r - To )^ + z^ ) ^ (3,41) 


The contributions from these stress terms are introduced into the 
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layer stress interpolation discussed in the previous section as follows; 
These terms are added to the 67 polynomial approximation to make it a 70 
term polynomial. To facilitate the programming of this condiotion, the $‘s 
from the top layer are introduced as shifting the jS's from the 

bottom layer to Hence, in Eqs. (3.39 - 3.4i), 

i = "71, j = 72 and k = 73 for the bottom layer 

while i = 16, j = 17 and k s 18 for the top layer 
ro = hole radius 
X X power of singularity 

Substituting Eq. (3.39) in equation Eq. (3.38), 

1 r 2^-1 

“^rz = j ^ ^ ''-''o ) + z ) * dr (3.42) 

ro . 

Solving this integral, the solution in terms of the Gauss 
Hypergeometric series is as follows 

Trz = - p i^i X 2 ( I ( ( s^ + 2 ^)^ - 2 ^) + rO z^ ■ •*• ( ( j® 

s^ + 2^ 


a a i 

^ s2 ' 

12 . (tt + 1) (0+1) j 

f 1 

7 ' 

^ . J2. i 

' 7 (7 + 1) 2 ! ' 

^ _2 2 ! 


(3 .431 


As this series converges quickly, only the first three terms of the 
series have been considered for the present formulation. 

At the hole boundary, 

CTr = 0 < as the boundary condition. (3.44) 

Substituting Eq. (3.43) and Eq. (3.44) in Eq.(3.37), 
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* S^ ( S^ + » ^ 5 gfl (o + i) (fl + 1) 5 7 7 4 

^ ^ rTTi) — s Cs + 2^) * 



£^ + 2^ 



...))) 


(3.45) 


Here, 


s = r-Po* 


a = 


1 

2 




Substituting Tpj from Eq.(3.40) into Eq.(3.38) 


o-z 


- i 6; 


+ z 


2f"l' 


- ;3. X 


4- Z 


2^1 


(3.46) 


where , ^ is the hypergeometric series 

For 2Fi^ ,a = |, 0 = 2±^and7 = | 

From Eq.(3.37) 


o-Q = r I z 0^ < ( r - ro )^ + ) ^ 


(3.47) 
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Finally; substituting Eq.(3.41) into Eq.(3.37) 

Trz = % J ^ )‘ 

S +2 

From Eq.<3.38) 

CTz = - Mk ((s^+z^)^ + + (s^+z^)^ ) 

(s2«2)^ 


(3.48) 


^ 1 00 < o+l)(g+D (,^2 2,2 , 2s^ 1 1 ^ 

2! 7 (7+i) > + I - 5 5 j j 




- h h ((-J— ; H. m 


1 s 


(s^+z^)* (s^+z^)* (s^+z^)* 


1 00 (a+iK0+l) ( 1 2s^ i s^ 

2! 7 (7+i) I „i‘ » ^l's ~ — 7 

(s^+z^)® 3(s^+z^)® (s^+z^)^ 


i)) 


(3.49) 


These stress cornponents which are in cylindrical co-ordinateS/ are 
transformed to the Cartesian system using the transformation matrix. 

The contribution from these singular terms is then introduced as 
the coefficients of - )S^2 layer i and the coefficient of - B^q 
for the layer i+l. The assembly procedure and the calculation of the H, G 
and the stiffness matrix K remains the same as has been discussed 
previously. 


33^ Determination of the Power of Singularity 

A finite element analysis has been carried out to determine the 
power of singularity (Yamada 1983 and Suresh 1987). For the axisymmetric 
case, the singular elements are a triangular sector spread across the hole 
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boundary as shown in the fig 3.4. In this approach, the singularity is 
determined by the use of isoparametric mapping. For the sake of 
completeness, this approach is explained in this section. The mapping is as 
follows .(fig .3 .5) 


and 





L 

2 


(3.50) 

(3.5i> 


e - + H2 ©2 + ^3 ^3 ( 3 . 52 ) 

where 

= I ( -7? + T?^) , H 2 = 1 - T?^ , H 3 = I ( T? + T?^) (3.53) 

are the shape functions and, s and 0 represent the polar coordinate radius 
and the angle of a generic point P, Sq being the radius of the common 
circular arc on which the nodes i, 2 and 3 are placed. 

A consistent isoparametric displacement field has been used as 
follows 

u - Us = [ (u^ - Us) + H2 (U2 - Us) + H 3 (U3 - Us) 3 ( ^ ) (3.54) 

where Us represents the displacement vector at the vertex S of the sector 
and Uj_(i=i, 2, 3) are the displacement vectors at the nodes on the circular 
arc of radius s = Sq. Substituting Eq. (3 .51) in Eq. (3.54) gives rise to a 
term in the expression for the displacement field. Hence, the singularity 
of stress and strain thus introduced is in the order of . 

The next step is to formulate the pertinent eigen value problem, 
which is done by solving the following axisymmetrio finite element prcAilem. 

The co-ordinates r and z, or the position of any point on the sector 


on the axisymmetrio plane is given by 
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r * a + s cose 

z = c sin 6 

as shown in the fig 3 . 5 . 


G.55) 

(3.56) 


The displacement vector is given by 



Expressing strains in terms of displacements, for the axisymmetric 

case 


€r 









8V2 

3z 

r 

it 

r 


I 

> 3r ; 32 ^ 


(3.58) 


The detailed derivation for the expressions for tp, and 7 pz 

are as given in reference SureshdSS?) which are obtained from Eqs. (3.50 - 
3.-57). This leads to 

c t } = CB3 Cup 


where CBl = ^ 5 ^ C X CBal + [8^3 3 ( 3 . 60 ) 

From the virtual work principle for the domain, the following 
expression can be had. 



({Xp Stp + a 2 3€-z + txg 3 fc 0 + Tp 2 STpj) 2 it r£jg dr dz 



32 


i 

= / 'fr a°r + ^ Tpz 3Cir 2i? r^n dO (3 .61) 

whS^e Peg is the radius of the centroid of the element and Pen is the 
radius of the node on the circular boundary with respect to the hole 
center. 

Substituting Eq.(3.60) into Eq.(3.6D; the characteristic equation in X 
is of the form 

[A] + X [B] + [C3 = 0 (3.62) 

which is then reduced to an eigen value problem for solution of X. 
The matrices LAD, CB] and CC] are given in reference SureshdSST). 

This eigen value problem has been solved for points along the hole 
boundary in order to determine the order of singularity as an angular 
variation. 
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CHAPTER 4 


RESULTS AND DISCUSSIONS 

This chapter presents the free-edge stresses in laminated composite 
plates subjected to a far field uniform strain £,y. The analysis has been 
carried out using hybrid finite element and singular hybrid element 
approach as discussed in the previous chapters. Results show the 

j;- 

ef fectivwness of the singular hybrid elements near the boundary. Two 
geometries have been considered for analysis: 

1. A symmetric crossply laminate with a central circular hole with 

[0/90]s and [90/03^ configurations. 

2. A symmetric crossply laminate with a central square hole with 

CO/9035 configuration. 

Details of the geometric and material data has been given in chapter IJ. 

4.1 Interlaminar Stresses Near the Circular Hole Boundary in Symmetric 
Composite Laminate 

As discussed in chapter-l there are situations where the composite 
laminates have holes, cut outs etc., which are a possible site for 
delamination. The inter-laminar stresses near the free edge of a circular 
hole in symmetric C0/90]s and [90/03s laminates have been investigated. 

Figure 4.1 shows a finite element discretization of the laminated 
composite plate with a central circular hole used in the present analysis. 
The laminate is subjected to a uniform far field strain in y^direction. As 
may be observed, only i/8^^ of the plate has been modelled due to the 
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symmetries of the plate about the x,y and z axes. The mesh is finer near 
the hole boundary and has a total of 85 elements. The elements are hybrid 
in nature with 67 fi's per ply (or 119 terms in all). As discussed earUer. 
only the elements near the hole boundary have been treated as singular 
hybrid elements with an additional three terms to model singularity. The 
singular hybrid finite element analysis has also been performed using a 
finer mesh having total of 170 elements as shown in figure 4.2. 

As a first step in formulating hybrid elements, the value of 
singularity of stresses near the hole boundary at the interface has been 
determined as explained in sec. 3.3, Fig 4.3 shows the variation of 
singularity values with the angle 6. As can be seen, the severity of 
singularity is the highest at Q = 0° and 9 = 90° as the material properties 
of the plies differ maximum at these two locations. At 9 = 45°, the reduced 
axisymmetric material properties do not differ and hence the singularity 
value is zero. These singularity values have been used in the construction 
of the stress terms in addition to the cubic variation, for the elements 
near the hole boundary. 

In what follows, these stresses have been normalised with a factor 
Sg as given by 

Sg=302 Uo/h Mpa 
Sg=43xi0^ Uo/h Psi 

Fig 4.4 - 4.6 show the interface normal stress, near the hole 
boundary obtained by using a coarse mesh of 85 elements. The variation is 
shown along radial lines at 6 = 0°, 6 = 45^ and 6 = 9D^. In the same figures, 
curves are also shown as obtained by Chahande(i989) without incorporating 
the singular elements and those of Raju and Crews(i982) obtained by 
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displacement formulation employing a very fine mesh(1900D degrees of 
freedom). As can be observed, the results with singular elements are much 
closer to those obtained by Raju and Crews(1982). The interface normal 
stress Cj will become negligible far' sway from the hole. However, in the 
present analysis, the rate of decay is a little less. This can be attributed 
to the coarse mesh away from the hole. These graphs prove beyond doubt 
the fact that the inclusion of the singular terms for the elements adjacent 
to the hole boundary improves the accuracy markedly and the value of the 
singularities in a simple axisymmetric analysis may be effectively used. 

Figs .(4 .7 - 4.9) show the plots of interface normal str'ess along the 
radial lines at 0 = 0°, e = 45° and 9 = 9D°. They show th« results obtained 
for a fine mesh with 170 elements (fig. 4.2) without singular elements and 
with singular elements adjacent to the hole boundary. It can be seen that 
the results are exactly same in the far field region but are different in 
terms of magnitude near the hole boundary. In the case of variation along 0 
= 0° line the present analysis shows very negligible interlaminar stress 
c 72 as is expected, However, at the free edge, there is a complete reversal 
in the nature of the stress. Comparing with Raju and Crews results, the 
compressive nature as determined by the analysis with singular elements is 
correct. Along 6 = 45° the stresses are compressive near the free edge 
and for 6 = 90°, they are tensile. It is worth noting that this is a possible 
location for delamination initiation and growth. 

Fig .4.10 shows the variation of the transverse Tpz shear stress in 
the radial direction along 0 = 0°, 0 = 45° and 0 = 90°. These stresses are 
zero at the hole boundary all in the three directions. This is due to the 
correct enforcement of the Tp^ = 0 boundary condition dLfl^ing the 
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derivauon of the einsuUr terms which is not eooarent in the results by 
Raju and Crews 

Figures 4.11 shows the hooo stress Oj variation near the hole 
boundary in the thickness direction at e = 0°, 6 = 45° and 6 = 90° One 

feature worth noui.-D is that this stress shows a distinct dicontinuity at 
the interface 

From the az through thickness distribution shown (fig 4.12), the 
following can be observed. The interlaminar stresses ctj have a large slope 
discontinuity at the interface as compared to the discontinuity observed 
by Raju and Crews This may be attributed to the fact that the interlayer 
traction continuity has been enforced in the present formulation. 

Fig 4 13 and Fig 4,14 show the variation of ctj and Tpz respectively 
in the three radial directions for the crossply laminate of the CSO/Ols 
configuration Except for the 0° direction, these stresses have changed in 
nature as compared to the C0/90]s case. This shows that the possibility of 
delaranation is also reversed. 

Figs 4 15 and 4 16 show the through thickness variations of the hoop 
stress Cq and Oz for the OO/Qls configuration. As in case of the [0/90]s 
results, the Oq distribution shows a discontinuity at the interface, while 
the Uz distribution shows a smooth variation with large slope. 

42 Edge Stresses Near the Sopuare Hole Boundary in Symmetric Composite 
Laminate. 

Similar to analysis of edge stresses near a circular hole in a 
symmetric laminate, singular hybrid finite element method analysis is 
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carried out in the case of a C0/OT3s laminate with a square hole. The 
laminate geometry and material properties are the same as in the previous 
case. The finite element discretization of the laminate is shown in figure 
4.17. The mesh is made finer near the hole and has a total number of 85 
elements. Analysis was also carried out using a fine mesh having 170 
elements (fig. 4. 18). However, there is one difference between the employment 
of singular terms in the construction of the hybrid elements adjacent to 
the hole. As the hole boundaries are straight (parallel to the x or y axes), 
the singularity is taken to be constant and equal to the value at either 
or 90^. 

The radial variation of the interlaminar stress cr^ and interlaminar 
shear stress Trz shown in figs .4. 19 - 4.20. These stresses have been 
plotted along the © = 0°, ©=45° and © = 90°. Here the normal stress cr^ at 
©=0° is negligible as was in the case of the circular hole problem, whereas 
it has significant positive value near S = 45 ° in contrast to the previous 
problem. At © = 90°, it follows the same trend as in the previous case and 
has a significant- positive value. These values seem to be much more than 
those in case of circular hole. However overall variation is the same. Thus 
square hole seems to have more tendency for delamination. 

The through thickness variation of the hoop stress and 

interlaminar normal stress are shown in figs .4 21 - 4.22. The hoop 
stress 170 as in the previous case, shows a distinct discontinuity at the 
interface but with a very high magnitude of the stress. Th interlaminar 
stress <7 zj due to the reasons discussed earlier, does not show 
discontinuity at the interface. 


It can be further concluded that the presense of significant 



44 


interlaminar normal stresses and interlaminar shesr stresses near the 
free edge, may act in conjunction, leading to mixed modes of failures. 
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nG.4-.t FEM DISCRETIZATION OF LAMINATE PLATE WITH CIRCULAR HOU 
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NO. OF ELEMENTS ■ 170 
NO. OF NODES ■ 565 

NO. OF TOTAL DEGREES 
OF FREEDOM « 50^ 

R ■ 5“ . b - 30" ; 1 - 6' 
THICKNESS CPLY) ■ 0.25“ 



FIG.4.2. FEM DISCRETIZATION OF LAMINATE PLATE 


WITH CIRCULAR H 
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Fig. 4. 4 radial distribution for 0/90 

with circular hole (85 elements) — e 


laminate 
= 0 ° 
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Fig. 4. 5 (Jz radial distribution for 0/90 
with circular hole (85 elements) - e 


laminate 
= 45'’ 


Fig. 4. 6 aj radial distribution for 0/90 laminate 
with circular hole (85 elements) - e = 90° 
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, through thickness distribution for 0/90 laminate 
with circular hole (170 elements) 
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Fig.4.12 Jz through thickness distribution for 0/90 laminate 
with a circular hole (170 elements) 


Fig. 4. 13 Gj radial distribution for 90/0 laminate 
with circular hole (85 elements) 
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Fig. 4. 15 a. 


through thickness distribution for 90/0 
with circular hole (85 elements) 


laminate 




Fig.4.16 <Ti through thickness distribution for 90/0 laminate 
with circular hole (85 elements) 





FIG.4.17. FEM DISCRETIZATION OF LAMINATE PLATE WITH SQUARE HOLE 
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FI6.4.18. FEM discretization OF LAMINATE PLATE WITH SQUARE HOLE 
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Fig. 4. 22 Gj through thickness distribution for 0/90 
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CHAPTER 5 


CONCLUSICINS AND SCOPE FOR FURTHER WORK 


5_1 Conclusions 

The present wcw^k deals with the edge effects in symmetric laminated 
composite plates with a central circular/square hole and subjected to a 
farfield uniaxial tensile load. Based on the results obtained in the 
previous chapter, the following conclusions may be arrived at : 

1. The singular hybrid finite element method shows a distinct improvement 
over the nonsingular hybrid method. This tool can be effectively used to 
study the behaviour of layered composites. 

2. For the CO/^ls configuration, the possibility of delamination exists in 
the direction perpendicular to the direction of loading, denoted as 0 = 90° 
in the previous discussion, where the transverse normal stress, cf^ has a 
high tensile value. The presense of high transverse shear stresses, viz. 
Tpz, in this direction may lead to mixed mode failures. 

3. In case of the square hole problem, the magnitudes of these stresses 
was much higher as compared to the circular hole case, indicating the 
higher tendancy for delamination. 

4. In all the cases, the interlaminar stresses have high magnitudes near the 
hole boundary and they die down quicldy, within a distance of approximately 
four times the laminate thickness away from the boundary. 

3J2 Scope For Future Work 


The present work, though satisfactory, has a few shortcomings whichj 



if overcome, may improve the results by a large extent while using lesser 
amount of computing time. 

1. Instead of the assumption of an axisymmetric case, a complete three- 
dimensional analysis should be attempted. The angular variation of material 
properties should be modelled in a better manner during the calculation of 
the power of singularity. 

2. The accurate enforcement of the traction free boundary conditions at 
the free edge should be made. 

3. New shape functions should be tried to achieve exact geometric mapping 
of the odd shaped elements near the hole boundary. 
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AP^ .^ 4 DIX - A 


The 69 i 3 stress assumption is given below : 

“ ^17 ^ 2^ ” ^ 54 ^ ” ^27"^^ ^28 y ^ 4 ^^2 * “ 5 ^ 29 ^ ^ 

+ ^ ^ 1^3- iSss) - 2 iSao^y + 1813 y^ + 1 1 (184 - 10 £ 7 > - 5 032 ^ 


+ c i ( 05 - - ^33] x^y + [ i ( 5 s - 053) - 3 03*3 xy^ + 0^9 y^ 

+ T ^ 03S + 033 X + 037 y + $33 X^ + $39 xy + 0^0 y^ 


+ ^ -3 < ^fis - ^8i5 + + 058 + 2 0^2 + 2 $35 ) x^y + 03i xy^ J 

(T y = $2^ + 0£2 X + t g( $7 - 0£j^) * ^ ^ ^23 ^ ^ 5 ^8 * ^€j} ~ 2 028^ 

+ t j( 09 ■ 052) “ g ^29^ y ■*■ ^24 ^ '*' ^ 5 ^^10 “ ^63^ ~ ^ ^31^ * ^ 


+ C 3 ( 0ii - 854) - 032 ^ xy + C I ( 0j^2 “ 855) “ 5 033 ^ y 

+ T 0^2 + 043 X + 044 y + 045 X + 0 ^ xy + 0 ^^ y + $^3 X y 


+ C -i < ^67 - 8i4 + 2 04 + 2 057 + 0i, + 054 ) Xy^ ) 
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f'xtf = & 2 S + <325 X + <327 y + & 2 e + <329 + <330 + <33^ x^ + 1832 

<333 Xy + ^34 y ^ ^49 '*’ < 3 s 0 ^ < 3 si y <3 s 2 ^ j ( 0gg 

~ <3j. 3 ^2 '*’ ^55 "*’ ^9 '*’ ^42 ^ ^ <338 '*’ ^37 ^ ^ Xy + 0gg y J- 

CT'yj = ti { I (i - T) (07 + 0g X + 0io + 0i2 y^ > + I ( ^ ^ ^ - 3 >^9 y 

+ J <1 - T) (-1 - 3 t) xy + I (1 + T) ( 05 Q + 05 i X + 052 ■*■ ^6S ^ 

+ g(l + T)(l + 3T) 052 y+g(l + T)(-i + 3T) 054 Xy 

+ i ( i - ) { ( 044 + 050 ) + ( 045 + 2 052 ) X + C -I ( - 0^3+ 02 

+0S£^ ■*■^< 347 * ^38^ 3 y + 043 X + C -g ^087 ~ 0i4^ - 3 C 04 + 0^7 ) 3 xy ^ ^ 
^'xx = { i ^ 1 + 08 y + 04 x^ + 05 y^ ) + I ( 1 - T ) ( 1 - 3 T ) 02 y 

+ J (i - t) (-1 - 3 T) 05 xy + I (1 + t) ( 054 + 055 X + 05 ^ X^ + 05 g y^ ) 

+ |(i + T)(i + 3T) 055 X + i(i + T)(-i + 3T) 05g Xy 

+ I ( i - ) I ( 035 + 05j^ ) + [ -| ( 055 - 0^3 + 09 + 052 ) + ( 038 



- ^ + <^39 + 2 /3s 3 ) y + C -f ( /3 c 8 - ^is > - 3 ( + 0«6 ) ] xy 

- «« } } 

ry 

Ct', = ^ I (1 - T)^ ( 2 + T ) ( 013 + 014 X + 015 y ) + i ( 1 + T ) ( 2 - T ) 

< <366+ <357 X +058 y) - i (i - T)^ (1 + T ) ( ( 02 + 09 ) + < 2 04 + 0ii > X 
+ ( 05 +2 012) y + J (i - T) (1 + T )^ [ ( 0££ + 052 ) + t 2 0£7 + 054 ) X 
+ ^ ®s® + 2 0££ ) y 3 



APPENDIX - B 


The standard bi-quadratic serendipity shape functions used for 8- 
noded quadrilateral elements are : 






N3(?,?7) = (i + ?)a-r?)(-i + ?- r?)/4 


N4(f,r?) = (i + S)(l-7?^) / 2 


Ns(S,y7) = (l + J)(l + r?)(-l + S-i-77)/4 


Ng(S,77) = (i-?^)(i + 7?>/2 


N7(S,77) = (i-S>(l + r?)(-i-? + r7)/4 


Ng(C,T?) = (l-S)(l + T?^>/2 



